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Defects in superfluid "^He, liigli-Tc superconductors, QCD colour superfluids and cosmic vortons 
can possess (anti)ferromagnetic cores, and their generalisations. In each case there is a second order 
parameter whose value is zero in the bulk which does not vanish in the core. We examine the pro- 
duction of defects in the simplest 1+1 dimensional scalar theory in which a second order parameter 
can take non-zero values in a defect core. We study in detail the effects of core condensation on the 
defect production mechanism. 

PACS numbers: 03.70.-|-k, 05.70.Fh, 03.65.Yz 



cn 

(N 
> 
O 
cn 



o . 

\o • 

o : 

Oh! 
I 

Oh, 



X 



I. INTRODUCTION 

Since phase transitions take place in a finite time, 
causality guarantees that correlation lengths remain fi- 
nite. If the symmetry breaking permits non-trivial homo- 
topy groups the frustration of the order parameter fields 
is resolved by the creation of topological defects to me- 
diate between the different ground states. Since defects 
are, in principle, observable, they provide an excellent 
experimental tool for determining how transitions occur. 

ZurekP, Q and Kibble^] originally suggested that 
causality alone is sufficient to bound the initial density 
of defects arising in a continuous transition, whether in 
condensed matter or quantum field theory. The analysis 
is very general, and depends on the fact that there is a 
maximum speed (e.g. the speed of light or the speed of 
sound) at which the system can become ordered. The 
ZK causal arguments can be quantified in many vari- 
ants, and we will not rehearse them here. It is sufficient 
to consider a system with critical temperature Tc, cooled 
through that temperature at a rate 



1 dT 



1 



(1) 



thereby defining the quench time tq. The prediction is 
that, if f is the defect separation at the time of defect 
formation, then 



(2) 



where tq is the relaxation time for short wavelength 
modes and ^o, also determined from the microscopic de- 
tails of the system, characterises the size of cold defects. 
The coefficient / is an undetermined efficiency factor, 
taken to be 0(1), but anticipated to be greater than 
unity. 



We term a the Zurek-Kibble (ZK) characteristic index. 
Its mean- field values are typically (7 = 1/3 for relativistic 
systems and (7=1/4 for strongly damped non- relativistic 
systems. Experiments on a range of condensed matter 
systems (superfluid ^He 0,0, high- Tc superconductors 
0,13], Josephson tunnel junctions (JTJs) gi'^^ 
suits that are commensurate with (|2Jl. 

This may seem paradoxical in that subsequent ana- 
lyti c approximations |lOl ITlj and numerical simulations 
[o. O ITsj have shown that, for simple systems, the 
scaling behaviour of ^ is understood, not so much in 
terms of causal bounds, but in terms of the instabilities 
of the time-dependent Ginzburg-Landau (TDGL) theory, 
whose dissipative behaviour controls which regime we are 
in. However, where scaling is appropriate it is found 
that causality arguments identify the correct engineering 
dimensions for the scaling behaviour of defect densities 
[H El and survives P|. 

For symmetry breaking of local gauge theories the re- 
sult ||2Jl is not complete because of the freezing in of long 
wavelength modes of the gauge fields 17]. However, for 
the high- Tc experiments of |6|, l^] and the JTJs of ^] ) 
the effect of this additional mechanism is small and we 
shall ignore it, given that it does not occur for the simple 
model that we shall solve numerically later. 

For each of the condensed matter systems listed above 
the relevant topological defects are vortex lines (strings) . 
Again, in the context of the early universe, strings (cos- 
mic strings) are the almost inevitable consequence of 
symmetry breaking in the most obvious extensions of the 
standard model for electroweak unification . Observa- 
tion has suggested the possible existence of cosmic strings 
J9j 20,, 21] but, as yet, the evidence is not compelling. 

Most simply, for both condensed matter and the early 
universe, the most studied and best understood strings 
arise from the breaking of a global or local U{1) symme- 
try, in which the order parameter can be represented by 
a single complex scalar field. In such cases the structure 
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of the vortex is quite simple, with a trivial core, in the 
centre of which the order parameter vanishes, restoring 
the U{1) symmetry there. That is, a string is a simple 
tube of false vacuum or ground-state, trapping flux if the 
symmetry is local. Although this is appropriate for the 
phases of the condensed matter systems for which exper- 
iments have been performed, it is a simplification since 
both superfiuid ^He and high-Tc superconductors have 
order parameters with several components. The unsta- 
ble ground state ('false vacuum') corresponds to $ = 0. 
We can separate the components of $ into two types 
<i> — {(p, ff) so that, in bulk, the (f) fields condense (i.e. are 
non-zero), whereas the fj fields do not. For any defect 
this characterises the situation in its exterior and, for a 
'normal' defect, both <j) and fj are zero in the core. How- 
ever, in some parts of parameter space it may require 
less energy to produce defects with 'abnormal' cores, in 
which 0, but the fj ^ 0. The existence of defects 
inside which the fj fields condense is not a question of 
topology, merely one of energetics. 

This is the case for superfluid ^He — B with its nine 
complex order parameters, for which there is experimen- 
tal evidence for superfluid '^He — B vortices with ferro- 



magnetic ife — A cores[23 123|- Also, this phenomenon 
has been predicted to occur in the idealised SO (5) model 
|24i l2fil | of high-temperature superconductivity, where, 
for appropriate doping, the cores of the conventional 
Abrikosov vortices are antiferromagnetic. There has been 
recent experimental evidence to support this theoretical 
picture^ IHIll 0, ^3 • 

This is not just a phenomenon of condensed matter 
physics. It has been argued that such strings 

also occur in the so-called colour superconducting phase 
of QCD, that is believed to be realized when the baryon 
density is a few times larger than nuclear density |34| . 
and hence in neutron stars. In fact, that such a phe- 
nomenon can occur in relativistic systems was first dis- 
cussed by Witten on looking at generalisations of 
cosmic strings. His model is the simplest of all that is 
compatible with early universe cosmology: consider a 
two component system described by two complex scalar 
fields (0, ?y) with an approximate U{2) symmetry. If the 
U (2) symmetry between fields <j) and 77 is explicitly bro- 
ken down to U{1) X ?7(1), the condensation might be 
energetically more favorable than 77 condensation, and 
the ground state will be given by {(f) 7^ and (77) — 0. 
As intimated above, this system permits (j) vortices char- 
acterized by the phase of (j) field varying by an integer 
multiple of 27r as one traverses a contour around the vor- 
tex core. Witten showed that, if the approximate U{2) 
symmetry is broken only weakly, the rj field may condense 
inside the core of the <j) string, breaking the corresponding 
C/(l) symmetry in the core. 

One consequence of this 'core-condensation' is that 
HE HE HE US it provides a way to stabilize a string loop 
against shrinking, as would normally be the case because 
of the string tension. The condensed matter counter- 



part to Witten's model of approximate C/(2) symmetry 
breaking is in two-component BECs |42j , for which both 
core condensation and this stabilisation may already have 
been seen. 

Most of the work cited above has addressed the static 
properties of defects with core condensation. In this pa- 
per we are interested in the dynamics of the creation of 
such defects. Insofar as the scaling behaviour of ^ is 
truly a consequence of causality we would expect it to 
be equally true for defects with and without core con- 
densation. However, it may not be so simple in that one 
consequence of core condensation is that the interactions 
between vortices can be drastically altered by the pres- 
ence of non-trivial cores . In the following sections we 
shall consider the simplest model permitting core conden- 
sation, to check whether the formation accords with the 
simple Kibble- Zurek scaling laws in tq of ((2Jl. 



II. THE 0(2) -» Z2 X ^2 ^2 MODEL IN ID 

Our model is a simplified version of more realistic dis- 
sipative systems, such as the Ginzburg-Landau descrip- 
tion of high-Tc superconductors. However, with the 
early universe in mind we extend it to a relativistic un- 
derdamped model when appropriate. 

Specifically, let us consider the symmetry breaking 

0(2) Z2X Z2 — > Z2 in one dimension, where ESB 
denotes explicit symmetry breaking by the introduction 
of terms in the action and SSB denotes spontaneous 
symmetry breaking. We consider a free energy of the 
form 



F 



dx 



[{d,(l>f + {d,7^f]+V{cjy,ij) 



(3) 



where 



y(0,,y) = ]^a{T)<j^^ + \a{TW + -^{<j^^+r^^f + '^^. (4) 

for which we have chosen to implement the 'phase transi- 
tion' explicitly through the Landau form a(T) = — a'(l — 
T/Tc), where a' > 0. The fact that there is no true tran- 
sition in one dimension is irrelevant for the rapid changes 
in temperature that we shall consider [T^ . 

When /3 = 1 F possesses 0{2) symmetry, and we can 
take (f) and rj to be the real and imaginary parts of a 
complex field $. There is a minor complication in that, 
if we fix the phase of $ to be zero, say, at infinity, making 
space , the winding number of the phase is a conserved 
charge. However, the configurations with non-zero charge 
are Skyrmions, and not the defects of relevance to the ZK 
scenario, of which none exist in the absence of topological 
charge when there are no fixed boundary conditions. 

The 0(2) is broken exphcitly to Z2 x Z2 when /? 7^ 1, 
which will be the case of interest for defect formation. In 
ID, the defects are kinks. For example, to convert F of 
© into an S0{5) model for high- Tc along the lines of ^3 
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we elevate 4> into the two-component complex field that 
couples to the electromagnetic field, and elevate 77 to the 
three-component Neel vector of antiferromagnetism. A 
similar transmutation would recreate a model like Wit- 
ten's. In this regard, we note that the potentials pre- 
sented in the literature are not quite identical in the way 
the different sectors couple. They only become so on im- 
posing non-linear constraints (e.g. cjr + ff = /^), as in 
[4^. At an effective level this can make little difference 
to the statics of defects (e.g. but for the dynamics 

we revert to the more fundamental linear sigma model of 
©, modelled on ^^'^ 01- Phase-transitions leading 
to kinks and walls with core condensation and displaying 
repulsive interactions were considered in the context of a 
S5X Z2 ^ S3X S2 theory in . In that case the authors 
were concerned mostly with the long-time behaviour of 
the system after the transition. In particular, significant 
changes to the scaling beaviour of domain wall networks 
were observed as well as the formation of defect lattices 
as the final product of the evolution. Here we will concen- 
trate on the period straight after the transition, namely 
on the process of defect formation, rather than on the 
long-time dynamics. 

The properties of the translation invariant configura- 
tions that minimize the free-energy above depend on the 
value of p. Defining (f>Q — —a/b and 7/q = —(3a/b we have 
that, for /3 > 1, 



F[0o,O] >F[0,r;o], 13 > I 

corresponding to an ground state with 4>{x) = 0. 
For < (3 <1 the opposite relation holds 

F[0o,O] <F[0,r;o], 0</3<l 



(5) 



and the preferred ground state configuration is now 
ri{x) = 0, (f>'^{x) = (/>Q. For /3 < 0, the case correspond- 
ing to Eq. |SJl holds again, as /3 is now negative and rj is 
therefore constrained to be zero. 




FIG. 1: Profiles of (j>{x) and rj{x) for several values of /3 > 0.5 
and a — —1.0 and 6 = 1.0 



turns out that, for this particular system, we were able 
to find such solutions analytically. 

Any static solution obeys the time-independent equa- 
tions of motion derived from the free-energy Eq. 



(9) 
(10) 



We start by assuming that 4>{x) has the usual kink-profile 



{x) — (/)o tanh(j:/^). 



(11) 



(6) where Z, the size of the core can now depend on (3. The 



trivial kink profile with ~ y^~2/a and 77 = is obvi- 
ously a solution for all values of f3. In order to find other 
possible solutions, we assume ri{x) ^ and replace the 
general kink profile (with arbitrary core size /) into the 
equation for 0, Eq. Solving in terms of rj we obtain: 



A. Kink solutions 

We will concentrate on the < /? < 1 regime and 
consider kink-like solutions that minimize the action in 
Eq. (O for the set of boundary conditions (f){—oo) = 
—4>o,ri{—oo) = and 4'i+^) = 4>0:Vi+'^) — 0- For 
any value of f3 in this region the free-energy is always 
stationary for the straightforward kink solution obtained 
by setting 



(x) = 00 tanh(Zfc2;), r]{x) 



0. 



(7) 



where Ik = (— is the kink size. The free-energy of 
this configuration is f3 independent and given by 



Fk 



(8) 



It is possible that for certain values of /3 other stable 
configurations exist with a lower free-energy ^3 ■ It 



r]{x) 00 1 + 



aP 



1/2 



cosh(a;/Z) ' 



(12) 



Demanding that this result satisfies the second equation, 
Eq. (|10|l we obtain an explicit expression for the kink 
core size 



I 



1 



1 



The full solution is then given by 
4>{x) — 00 tanh(a;/Z) 
f]{x) = ±00^2/3-1- 



1 



cosh{x/l) ' 



(13) 

(14) 
(15) 



Note that we have included the two possible signs for 
the 7] field. In general, for each choice of profile, there 
will be two allowed symmetric configurations for rj. We 
will refer to solutions for which interpolates between 
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FIG. 2: The free energy of the most stable kink configuration FIG. 3: Maximum value of r;(x) and core-length of the most 
against /3, normalized by the energy of the rj = kink. stable configuration versus f3. 



—00 £ind (j)o as x goes from — oo to oo as kinks. These 
can be positively of negatively charged, depending on 
whether rj is positive or negative. Conversely, anti-kinks 
will have (/>(— cx)) = and 0(oo) = — 0o- We will say 
an anti-kink has negative charge if 77 > and positive 
charge otherwise. This convention for the charge signs for 
kinks and anti-kinks corresponds to the winding direction 
of the solution in the (</>,??) plane, mimicking the U{1) 
charge definition in the /3 — > 1 limit. As we will see 
below, the interaction potential between kinks and anti- 
kinks will depend on their charges, as defined above. 

The free-energy of the non-trivial kink profile can be 
evaluated explicitly, being given by 



F^FkV2-2(3i^-+pj , (16) 

where Fk the energy of the 77 = solution as defined 
in Eq. ©. Such a kink is rather like a a slice through 
a configuration known as a 'dark-bright vector soliton' 
in two-component BEC which consists of a domain 
wall formed by one condensate with the second conden- 
sate confined to the wall's centre. The difference lies in 
the fact that the BEC condensate obeys non-dissipative 
Gross-Pitaevskii equations, rather than the (relativis- 
tic) dissipative TDGL equation appropriate to condensed 
matter and the early universe. 

The above expressions for the fields' profiles and the 
free-energy of the non-trivial solution are defined only 
for f3 > (3c, where /3c = 0.5, a value independent of the 
parameters of the theory. For /3 above the critical value 
we have that F < Fk, the free-energy of the non-trivial 
solution is smaller than that of the 77 = kink. In this re- 
gion of parameter space the configuration in Eqs. (|11I12|) 
is energetically more stable and should be more likely to 
form as the outcome of a phase-transition. As (3 ^ /3+ 
this solution converges smoothly to the trivial kink pro- 
file, with Z ^ ^fc, 77 ^ and F -> Ffe. For /? < Pk, Eq. Q) 



corresponds to the single valid solution. These points are 
illustrated in Figs. In Fig. ^ we show several field 
profiles in the (3 > (3c regime. As (3 approaches its critical 
value the size of the kink core decreases tending to Ik , and 
the magnitude of the peak of 77 goes to zero. In Fig.|21we 
have the free energy of the most stable configuration, for 
(3 in the interval (0, 1). Finally, in Fig. Owe show, for the 
most stable configuration, the maximum value of rj and 
the size of the kink core. As expected, both these and 
the free-energy are constant below (3c- It is interesting 
to note that while all these quantities are continuous at 
(3c, their derivatives behave in different ways. In partic- 
ular, F'{P) is continuous throughout the interval (0,1), 
whereas the core size has finite discontinuous derivative 
at (3c and 77max has diverging derivative at 



B. Kink/anti-kink interaction 

The problem we want to address here is to determine 
what happens when a kink and an anti-kink are put in 
the vicinity of each other. For the trivial profiles corre- 
sponding to (3 < 0.5, it is well known that a kink and 
anti-kink attract each other with a force that decreases 
exponentially with the separation of the cores. In a set- 
ting where kinks and anti-kinks form as a consequence of 
a phase-transition, this allows for pair annihilation pro- 
cesses. In our case, because of the presence of the extra 
field, kinks can have two distinct charges and it is rea- 
sonable to expect these to play a role in their interaction. 
If we imagine a kink/anti-kink pair with the same charge 
in the limit /3 ~ 1, the result is a configuration close to a 
"springy" i7(l) winding. The U{1) spring will tend nat- 
urally to "stretch" itself, corresponding, in terms of the 
kink/anti-kink pair, to a repulsive interaction. The op- 
posite reasoning suggests that kinks and anti-kinks with 
symmetric charges should attract each other. 
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r/l(p) r/l(p) 



FIG. 4: The gradient flow interaction potential between a kink 
and an anti-kink with opposite charges, for several choices of 
/3. The kink/anti-kink separation is given in units of core size. 



FIG. 5: Interaction potential for a kink/anti-kink pair with 
identical charges, for a selection of values of As in Fig. |1| 
the distance is in terms of the static solution core size. 



In order to check this scenario and obtain effective 
interaction potentials we followed the gradient flow ap- 
proximation |47j| . This consists in setting up an initial 
kink/anti-kink configuration and evolving it numerically 
according to the first-order equation of motion corre- 
sponding to the free energy in Eq. 

dV 

dt<Pa = dl^a - TT-, (17) 
OCpa 

where the index a runs from 1 to 2 with (pi — (f> and 
(j)2 = rj- This forces the solution to the be locally on 
a minimal energy configuration for all t, in the sense 
that its trajectory should follow the "valleys" of the free- 
energy landscape. At given simulation times we define 
the defect separation r as the distance between the two 
zeros of cj) and we measure the total free-energy of the 
system. A potential V{r) for this "moduli" is defined by 
subtracting the energy of two isolated kinks to the energy 
measured for each value of r in the simulation. 

In Figs, ^and^lwe show V{r) for several values of /3, 
for pairs with opposite and identical charges respectively. 

As expected, the force between opposite charged kinks 
and anti-kinks is always attractive. As /3 increases and 
the kinks become less localized, the steepness of the po- 
tential decreases indicating a weakening of the interac- 
tion force. A kink/anti-kink pair with identical charges 
displays, on the contrary, a repulsive interaction for dis- 
tances above a certain threshold. For smaller distances 
the interaction seems to be always attractive, though re- 
sults below the core size, r < l{j3), should not be trusted 
as the "moduli" becomes ill-defined. As (3 increases, the 
height of the potential maximum becomes larger and its 
location moves towards smaller values of r. Though for 
values of (3 just above the 1/2 threshold the size of the 
potential barrier is likely to have negligible effects, as 
/3 — > 1 the kinetic energy needed to surpass the barrier 



will become increasingly high. As a consequence of the 
repulsive nature of their interaction, same-charge kinks 
may be prevented from colliding and we can expect a 
regime where pair annihilation will be suppressed. As 
we will see, this will lead to qualitative changes in defect 
production in a typical phase transition. 

III. NUMERICAL SIMULATIONS OF DEFECT 
FORMATION 

We now proceed to numerically evaluate the forma- 
tion of these various types of defects by using a Langevin 
equation to simulate a series of quenched transitions. Our 
main goal will be to check whether the types of scaling 
predicted by the Kibble-Zurek scenario in Q are com- 
patible with the presence of non-trivial defects as de- 
scribed in the previous section. We will evolve the gen- 
eral second-order Langevin equation 

dV 

d^C^a - dl<i>a + O? dt<t>a + TJ— = "Ca, (18) 

where the index a runs from 1 to 2 with (pi = (f> and 
(p2 — rj. is a gaussian noise term obeying 

{Ca{x' ,t')a{x,t)) ^T5{x' ~x)5{t' ~t)5ab, (Ca(x,t)) =0. 

(19) 

and a measures both the amplitude of the noise and 
of the dissipative first-order time derivative term. The 
relation between these two terms ensures that the 
fluctuation-dissipation theorem is satisfied and guaran- 
tees that for very large times thermal equilibrium at tem- 
perature r/2 is reached. A simpler version of this model 
with one single field has been used successfull y in several 
studies of defect formation in 1-1-1 dimensions jl3.l48l |: we 
refer the reader to these for more a detailed discussion of 
the model and its numerical implementation. 
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Our goal here is to simulate a "phase-transition" with 
general quench rate Tq^ . We start at high-temperature 
with a > and, after allowing the fields to equilibrate, 
we decrease a(T) with T changing according to Eq. 
As T goes below Tc, a series of defects forms separating 
regions in alternating vacua. Finally, a(T) settles to a 
constant negative value and the defect population enters 
a period of relative slow evolution, with occasional pairs 
of defects and anti-defects annihilating. At this point, 
we measure the defect density p = 1/^ by counting the 
number of kinks and anti-kinks in the field configura- 
tion. The defects are identified as zero crossings of the (p 
field, suitably coarse-grained over a few lattice sites. The 
coarse-graining is particularly relevant for high values of 
f3, when the size of the defect core becomes considerably 
large, and small fluctuations around zero might lead to 
defect over-counting. 

By repeating the above process for several values of 
the quench rate we obtain p{tq) and fit the results to 
a power-law Atq"^ . Note that depending on the magni- 
tude of the parameter a, Eq. H18I) may describe either a 
relativistic system (small a) or one where dissipation is 
dominant (large a). The former type of setting will be 
typical of high-energy and cosmological systems, whereas 
the later will correspond to condensed-matter systems. 
As discussed in Section U we should obtain, according 
to the standard KZ predictions, cr = 1/4 for the over- 
damped regime and cr = 1/3 for the under-damped case. 

In Fig. El we show cr as a function of (3 for two choices 
of q;, representative of the two types of regime. In both 
cases we set a — 1, b = 1 and F = 0.01, and we vary the 
quench rate as tq = 2", n — 1,2,.., 9. The final defect 
densities are obtained by averaging over ten independent 
realizations. Production runs were performed on periodic 
lattices with 16000 points and space and time steps of 
6x = 0.125 and St — 0.1. We found that coarse-graining 
(j) over eight lattice points was enough to eliminate any 
unwanted fluctuations of the field. The error-bars, not 
shown in the plot for clarity, were obtained by calculating 
the standard deviation over the ensemble of 10 sample 
runs, and are of order of 10% of the values observed for 
a. 

Starting with the high-dissipation results with a = 1.0, 
we see that within error-bars, the scaling exponent is ap- 
proximately constant for /3 < 0.5. The value measured 
a ~ 0.28, coincides with that obtained for similar sim- 
ulation parameters in the case of a single field system 
|48l |. One should note that the method used to deter- 
mine the final defect density is known to lead to a slight 
over-estimate in cr when compared to other, computa- 
tionally more demanding, approaches (for a discussion 
see [12 )• This leads to the deviation from the theoreti- 
cal prediction cr ~ 1/4 for the over-damped regime. With 
this caveat in mind, and remembering that the transition 
between the under and over-damped regimes is continu- 
ous, we will proceed as in and use the value cr ^ 0.30 
to distinguish between the two types of behaviour. In 
the following, all values of cr should be taken in refer- 
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FIG. 6; The scaling exponent ct as a function of f5 for 0.4 < 
f3 < 1, in the under-damped (a — 0.6) and over-damped 
{a = 1.0) regimes. The error-bars, not shown for clarity, are 
of order of 10% of the values measured. 



ence to those obtained in the low-/3 parameter region. 
As is clear in Fig. El as /? increases, the scaling expo- 
nent goes down signaling a departure from the canonical 
KZ-scaling behaviour, moving deeper into the dissipative 
regime. For values of (3 in the vicinity of 0.80, cr reaches a 
stable plateau, remaining constant (within error-bars) up 
to P ^ 0.97 when a sharp increase takes place. The re- 
sults for the relativistic case a — 0.6 show a very similar 
pattern, with a steadily decreasing towards a flat regime, 
followed by a sudden rise as /3 ^ 1. 

The initial period of slow decrease in a for values of /3 
up to 0.7 — 0.8 can be related to the fact that the effec- 
tive physical consequences of the dissipation should be 
measured in terms of the relevant masses of the problem. 
The relaxation time scale that determines the type of be- 
haviour (dissipative vs relativistic) is given by the inverse 
of o? jm^ where m stands for the lowest mass in the the- 
ory. As /3 increases and T < Tc^ the mass of one of the 
excitations of the vacuum approaches zero. When 13—1 
and the original Lagrangian becomes explicitly U{1) sym- 
metric, this degree of freedom corresponds to the mass- 
less Goldstone boson. More specifically, the masses of the 
vacuum excitations for the potential in Eq. are given 
by ra\ — —2a and m| = — a(l — /3). Clearly, m2 takes 
increasingly small values as /3 approaches 1. This leads 
to a rise in the largest effective dissipation scale a'^/m2, 
shifting the system towards the over-damped limit and 
decreasing the value of cr. In order to check this we sim- 
ulated a one field relativistic system with a varying mass 
term. We observed as expected, that as the mass was 
reduced the value of cr decreased, with the system dis- 
playing dissipative effects. 

How is this trend towards lower values of the scaling 
exponent stopped? To understand the overall behaviour 
of cr as displayed is Fig. O one must keep in mind that as 
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FIG. 7: Fraction of same-charge kink/anti-kink pairs 
(crosses) and opposite charged pairs (squares) as a function of 
j3. For reference we also show the fraction of negative charge 
kinks (diamonds) . The bottom plot corresponds to the under- 
damped regime a = 0.6 and the top plot to the over-damped 
case Q = 1. 



f3 takes larger values, the interaction between kinks and 
anti-kinks changes considerably. As discussed in Section 
IliijI for large separations, there is a repulsive interac- 
tion between kinks and anti-kinks with the same charge. 
Note that the effect of the repulsive interaction should 
not be felt immediately after (3 > 0.5, as the height of 
the potential barrier is initially very small. Nevertheless, 
since the energy needed to overcome the repulsive barrier 
increases with (3, we expect kink/anti-kink annihilation 
to be inhibited for larger values of this parameter. This 
should lead to survival of a larger number of kink pairs, 
increasing the final number of defects. What would be 
the signature of this mechanism in the scaling of the fi- 
nal defect density with quench rate? Clearly the effect 
would be felt more strongly for fast quenches, where de- 
fect densities are higher and annihilation is likely to play 
a bigger role. An increase in kink survival rates for low 
values of tq should lead to a steeper distribution p{tq), 
that is, larger values of the exponent a. 

The mechanism above provides an explanation to why 
the decrease in a is halted for values of f3 roughly halfway 
between 0.5 and 1.0. In order to confirm this scenario 
we identified neighbouring pairs of kinks and anti-kinks 
in the simulation and determined their charges by look- 
ing at the sign of ij in their cores. Using the knowl- 
edge of the spatial distribution of charges at the final 
time of each simulation, we counted the number of de- 
fect pairs with equal and opposite charge respectively. 
If same-charge pairs are being prevented from annihilat- 
ing, then their numbers should be in excess of those of 
opposite charged pairs. In Fig. |7| we compare the total 
number (i.e. summed over all quenches) of both types of 
kink/anti-kink pairs as a function of in the dissipative 
and relativistic cases. The results match our expectations 



well, with the fraction of pairs with equal charge deviat- 
ing little from 0.5 up to values of ~ 0.7. Above that 
threshold there is a decrease in the percentage of pairs 
with attracting charges, signaling the survival of equal 
charge pairs due to the rise in the repulsive potential bar- 
rier. We also confirmed that this effect is more marked 
for fast quenches, as discussed above. For reference, the 
plots include the data of the fraction of pairs with nega- 
tive charge as well. As expected, as a consequence of the 
symmetry of the theory, this quantity remains equal to 
0.5 within error-bars throughout the whole range of p. 

Finally we will focus on the apparently anomalous be- 
haviour of the system for values of /3 very near 1, in which 
limit there is no topological charge. For the three points 
with highest /9, i.e. for /3 > 0.96, there is a marked in- 
crease in the exponent tr, accompanied by a decrease in 
the fraction of same-charge pairs. A closer look at the 
data reveals that the change in relative amount of types 
of kink/anti-kink pairs is caused exclusively by the num- 
ber of pairs with opposite charge going up by a consid- 
erable amount. We note that the results for these three 
values of (3 should be taken with care since we are in 
a regime where the kink core size becomes increasingly 
large. In particular, for fast quenches with high final 
defect densities, the defect-defect separation becomes of 
the order of magnitude of the core size. Taking the most 
extreme case (3 = 0.99, we find for the simulation param- 
eters / = 10, which implies that each defect has a spatial 
extension of the order of 20. Our simulation box should 
in principle not be able to accommodate more than 100 of 
these kinks, and for fast quenches the number observed 
is indeed very close (for example 80 in the relativistic 
case, for the smallest tq). In other words, the correla- 
tion length leading to domain formation at the transition 
is close to the core kink size. These cases should be un- 
derstood as being effectively in the U{1) regime, with 
winding springs being formed with no spatially localized 
core. The dynamics of these objects should be described 
in terms of unwinding processes and this might lead to a 
slower annihilation rate than traditional kink/anti-kink 
collisions. In terms of the interaction potentials in Fig.^ 
this corresponds to the flattening of the attraction po- 
tential as /3 ^ 1. A study of the time-scales involved 
in the departure from equilibrium in such case (using for 
example the techniques developed in jS^) could also shed 
some light on the dynamics of these systems. A detailed 
analysis of the unwinding process and its role in defect 
production in a phase-transition will be the focus of a fu- 
ture publication dedicated exclusively to the U{1) theory 
in 1 + 1 dimensions 15011. 



IV. CONCLUSIONS 

We studied kink formation in 1-1-1 dimensions, in a 
scalar theory where an added field can lead to conden- 
sation in the defect core. We obtained explicit expres- 
sions for the non-trivial defect profiles, and determined 
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the region of parameter space for which these are sta- 
ble. This introduces an extra degree of complexity in 
the theory as kinks acquire a charge defined in terms of 
the sign of the core condensate. As a consequence, the 
process of defect production at a non-equilibrium phase 
transition suffers a qualitative change, with annihilation 
of kinks and anti-kinks being suppressed due to the re- 
pulsive nature of their interaction. The degree to which 
this takes place depends on the parameters of the po- 
tential. As a consequence we observed deviations from 
traditional Kibble-Zurek scaling, with the value of the 
scaling parameter showing a clear, measurable signature 
of this effect. 

These results open interesting possibilities for the case 
of higher dimensional theories. Though there seems to 
be good evidence by now in favour of the stabilit y o f 
vortons in U{1) x U{1) theories in 3-fl dimensions 
the question of formation still remains unresolved. Vor- 
tons are very fragile objects making a numerical study 
of their production in a phase-transition impossible, the 



size of the domains required being severely constrained 
by computational limitations. The model described in 
this article can be generalized to 2+1 dimensions and an 
extra condensate field introduced. Such Z2 x U{1) the- 
ory would display domain walls formed by kinks of the 
Z2 field, with the U{1) field condensing in their cores and 
leading to vorton-like solutions. The lower space-time di- 
mensionality of this model could be make numerical sim- 
ulations of vorton formation possible - a line of research 
we will be following in a future publication. 
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